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refrigerant used in the system. The design of the system using PCS needs the quantitative information 
about the thermal and fluidic behaviors of PCS. Here we provide an overview of the characteristics of two 
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1. Introduction 
With the rapid development of the economy and human 
society, the energy and environmental issues have become more 
* Corresponding author. Tel.: +86 21 34205505; fax: +86 21 34206814. and more remarkable focuses because they will finally threat the 
E-mail address: zhangp@sjtu.edu.cn (P. Zhang). survival of the human being if we do not take any actions to deal 
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Nomenclature 


parameter in Eq. (3) 
parameter in Eq. (2) 

specific heat capacity (kJ/(kg K)) 
diameter (m) 

shear rate (s7!) 

Fanning friction factor 
volumetric flow rate (m?/s) 
heat transfer coefficient (W/(m? K)) 
Hedstrom number 

enthalpy change (kJ) 

length (m) 

latent heat (kJ/kg) 

k thermal conductivity (W/(m K)) 
K, K’ fluid consistency coefficient 

m mass (kg) 

M molecular weight (g/mol) 

n flow behavior index 

N hydration number 

Nu Nusselt number 

Pe particle Pelect number 

AP pressure drop (Pa) 

q heat flux (W/m?) 

Q heat (W) 

r radius (m) 

radius (m) 

Reynolds number 

the modified Reynolds number 
Stephan number 

temperature (K) 

temperature difference (K) 
velocity (m/s) 

mean velocity (m/s) 

mass fraction (wt%) 

power (W) 

mass fraction 

molar fraction 


i) 


oO 


-—~prrTaArAWA QDs 
an) 


x <2) 
x SPCR PIL RRs 


Greek letters 


a thermal diffusivity (m/s) 
y shear rate (s7!) 

À Darcy friction factor 

H viscosity (Pa s) 

Ha apparent viscosity (Pa s) 
Ho viscosity for Bingham equation (Pa s) 
p density (kg/m?) 

T shear stress (Pa) 

To yield stress (Pa) 

(o volume fraction 
Subscripts 

0 initial 

a apparent 

b bulk 

B Bingham fluid 

e effective 

f fluid 

H hydrate 


i inlet 

l liquid 

MR Merzner and Reed 
Newtonian 

m phase change material 
particles 
slurry 
wall 
distance from the inlet 


x~evrvyvaZ 


with them. The countermeasures to the energy issue may include 
improving the energy efficiency, making use of the renewable 
energy and waste heat, and so on. One of the heavy energy 
consumers is the refrigeration and air conditioning; meanwhile, 
the ozone layer depletion and greenhouse gases emission of the 
environmental issue pose big concerns to the air conditioners and 
refrigerators where a huge amount of refrigerant is used, among 
which CFCs and HCFCs are the primary ones. The phase-out of CFCs 
and HCFCs according to the Montreal Protocol urges the 
researchers to find out the environment-friendly new substitu- 
tions for refrigerants or new technology to reduce the quantity of 
the refrigerant used in the systems. One approach to reduce the 
quantity of the refrigerant used in the circuit is to use the 
secondary refrigeration loop, which is responsible for distributing 
the cooling power. The main merits of such system are the quantity 
of the refrigerant involved is significantly reduced and the energy 
efficiency can be possibly improved by using the proper working 
medium in the secondary loop although the additional circulating 
device-pump is necessary. 

The essence of the application of PCS to the secondary 
refrigeration loop is the introduction of the energy storage into 
the system, which releases the cooling power and enables that the 
primary refrigeration system does not need to work or works 
under a reduced load condition when the cooling power is needed, 
so that the electricity peak shaving can be readily realized. To 
increase the energy storage capacity, two-phase slurry fluid is 
preferred which is inspired by the fact that the latent heat storage 
associated with solid-liquid phase change is more efficient than 
the sensible heat storage of the single liquid phase. Apparently, 
selecting the appropriate working fluid for the secondary 
refrigeration loop plays a very critical role in such system. 

In general, there are principally four kinds of two-phase fluids 
which can potentially be used as the secondary refrigerants [1], 
which have been named as the latent functionally thermal fluids 
(LFTF) [2]. Among these fluids, although ice slurry has been 
subjected to very intensive research and can be a very straightfor- 
ward choice for the secondary refrigeration loop, the practical 
application is still limited due to the facts that there are still some 
weak points need to be solved, such as, the requirement for reliable 
ice slurry generator and ice concentration controller, and the 
requirement for improving the refrigerating efficiency of the ice 
slurry generation system because of significant subcooling neces- 
sary for ice nucleation. Recently, significant attention has been paid 
to the microencapsulated phase change slurry (MPCS) and semi- 
clathrate hydrate slurry (CHS), which have potential to be used inthe 
secondary refrigeration and air conditioning loop systems. The main 
merits of MPCS and CHS over ice slurry are as follows: (1) the phase 
change temperature range is well fitted for air conditioning system, 
i.e. 5-12 °C, by properly selecting the phase change material (PCM) 
or by varying the salt concentration of the aqueous solution in the 
case of CHS; while the phase change temperature for ice slurry is 
0°C, lower than needed for the room conditioning; (2) the slurry 
particle size of MPCS or CHS can be smaller than that of ice slurry, 
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Table 1 
Comparison of the typical materials used for PCS. 
Material Melting Density Latent 
temperature (kg/m?) heat (kJ/kg) 
CC) 
MPCM (C16H34) 18 835/776 237 
(solid/liquid) 
CHS (TBAB, 11.8 1080 193 
(C4Ho)4NBr-26H20) 
Ice 0 920 333 


which results in smaller frictional pressure loss for the same mass 
flow rate and less risk of clogging the transportation pipes; (3) the 
energy efficiency of the refrigerating system for MPCS and CHS is 
higher than that for ice slurry because of higher phase change 
temperature. The comparison of three typical materials for PCS is 
briefly shown in Table 1. Although the latent heat of ice associated 
with solid-liquid phase change is the largest, it is dimmed by the 
drawbacks stated above. 

This paper presents a literature review of the research progress 
on the two-phase secondary loop by employing MPCS and CHS as 
the working fluids. The review has been arranged into two big 
sections regarding MPCS and CHS, respectively. Each section 
begins with the selection of the proper materials for PCS and is 
followed by the determining of the thermal properties of the 
slurry; the fluid flow and heat transfer characteristics of PCS are 
further discussed as they are very essential for the system design 
and optimization; at the end of each section, some practical 
applications of PCS to the refrigeration and air conditioning are 
mainly introduced and some problems are also proposed for 
further research in the summary. 


2. Microencapsulated phase change material slurry (MPCS) 


2.1. The material and fabrication of microencapsulated phase change 
material (MPCM) 


Microencapsulation is a kind of tiny particles (about 0.1- 
1000 um in diameter) in which the core material is enveloped by a 
thin layer of shell. The core material can be solid, liquid, gas or even 
the mixture, and the shell material is generally natural or 
synthesized polymer of about several micrometers in thickness. 
This technology has been widely used in such fields as medicine, 
chemical engineering, etc. Recently, it has been extended to the 
preparation of the microencapsulated phase change material 
(MPCM), by which several inherent problems, such as low thermal 
conductivity, big volume change of the organic PCM and the 
thermal instability of the inorganic PCM, can be possibly overcome. 
When MPCM is dispersed into the carrier fluids, it forms a kind of 
two-phase suspension-MPCS and the application of MPCS to the 


oR 


energy system becomes very simple. Although the carrier fluid 
should have high thermal conductivity and large specific heat 
capacity, water is generally used as it has no obvious negative 
effect on the compatibility of MPCM and is easy to handle. 

There are several methods which can be used to produce 
microencapsulated particles, such as spray-drying and fluidized 
bed process for physical approaches; and there are also chemical 
processes which can be used for microencapsulation, such as 
simple and complex coacervation, in situ polymerization and so on. 
The in situ chemical process is often used in the laboratory since it 
can produce the microencapsulated particles with good quality. 
Although the diameter of MPCM depends on the fabrication 
method, the typical size is in the range of micrometers. Shown in 
Fig. 1 are the typical SEM pictures of MPCM [3]. It can be seen from 
the surface morphology that the shape of MPCM is uniformly 
spherical and the surface is very smooth without edges and dents, 
which implies that fluidity is good if the MPCM is dispersed into 
the carrier fluids, such as water. 

The intensive effort for developing MPCM can be traced back to 
two or three decades ago in America [4]. However, much attention 
was not paid to the application of MPCM to such fields as thermal 
system, building material, textile application and so on, until 
recently. 

There are many kinds of MPCMs which can potentially be used for 
the secondary loop purpose. However, before it can be used, one kind 
of MPCM should generally meet several requirements for both core 
and shell materials, which include the follows: (1) core material: 
suitable temperature range, large latent heat, reasonable thermal 
conductivity, low volume change, low reactivity, and so on; (2) shell 
material: good sealing tightness, endurance, good elastic strength, 
water and fire resistance, and so on. It is actually almost impossible 
to find a MPCM which can meet all the above requirements. 

The core material used for MPCM is commonly organic material 
considering the technical feasibility, and the shell material can be 
urea-formaldehyde resin, melamine-formaldehyde resin and 
even some other resins depending on the chemical process 
involved. It is obvious that the paraffin is the best candidate as 
the core material due to its stable chemical properties, suitable 
thermal properties and low cost. The thermal properties of the 
paraffin depend on the carbon atom number in the molecule, and 
the melting temperature ranges very wide from —57 °C for n- 
Octane (CsH18) to 80°C for n-Tetracontane (C4oHs2). The latent 
heat associated with the solid-liquid phase change of paraffin is 
quite large. Table 2 shows the technical data for n-alkanes which 
are the main contents of the paraffin. However, the heat transfer is 
always limited by its low thermal conductivity during heat 
retrieval. One solution is making the paraffin into microparticles by 
utilizing microencapsulation technology. 

Alkan et al. [5] have prepared MPCM with docosane as the core 
material and polymethylmethacrylate (PMMA) as the shell 


Fig. 1. Typical SEM pictures of MPCM made by different methods [3]. (a) By spray-drying and (b) by coacervation. 
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Table 2 


The technical data for n-alkanes. (The values of the thermal conductivity and density may vary with temperature). 


Carbon Number Name Formula (mol weight) Melting 


Latent Heat Thermal Conductivity Density (kg/m?) 


Temperature (°C) (kJ/kg) (W/(mK)) solid/liquid solid/liquid 

c(8) n-Octane CgHjg (114.23) —57 (—56.5) 181 —/0.128 —/[703.6 
c(10) n-Decane CjoH22 (142.28) —30 201 —/0.1349 —/[730 

c(12) n-Dodecane C12H2¢ (170.34) —9.5 (-12) 216 0.21/0.14 —/750 

c(14) n-Tetradecane C,4H39 (198.39) 6 226 0.21/0.15 825/771 
c(16) n-Hexadecane C1¢H34 (226.44) 18 237 0.21/0.151 835/776 
c(18) n-Octadecane CgH33 (254.50) 28 243 0.35/0.151 840/778 
C(20) n-Eicosane C20H42 (282.55) 37 247 0.23/0.151 856/780 
C(22) n-Docosane C22H46 (310.61) 44.5 249 0.22/0.151 864/780 
C(24) n-Tetracosane C24Hso (338.66) 52 (51) 255 0.25/0.15 897/774 


material by using emulsion polymerization. The characterizations 
of MPCM were also carried out by using SEM and FT-IR for the 
morphology and chemical contents. Differential Scanning Calori- 
metry (DSC) and Thermo-Gravimetric Analysis (TGA) were applied 
for determining the thermal properties and thermal stability. The 
results indicated that the MPCM was very suitable for thermal 
energy storage for space heating purpose since its phase change 
temperature was around 40 °C. 

A systematical investigation of the synthesis of MPCM with 
paraffin wax with high melting temperature as the core material and 
urea-—formaldehyde as the shell was carried out by Jin et al. [6]. They 
focused on the effects of the synthesis conditions on the morphology 
and thermal properties of MPCM and found that the particle size of 
MPCM decreased from 20 to 5 um as the amount of the emulsifier 
increased. The quantity of the shell material used in the synthesis 
was very crucial in deciding the toughness of the shell of MPCM. It 
was shown that the MPCM particles made under certain conditions 
might experience leakage when they were subjected to the thermal 
cycles with the temperature higher than the melting temperature. 
Consequently, the core material (paraffin wax) re-coalesced, which 
implied the failure of the MPCM particles. The pipe might be clogged 
if the leakage of MPCM occurs since the core material may stick 
together to form big lumps. However, this could be avoided by 
properly increasing the amount of the shell material used in 
synthesis process. DSC tests indicated that the melting temperature 
of MPCM was little affected and the latent heat decreased almost 
linearly with the core content. 

A simple method to make MPCM has been proposed by Ai et al. 
[7]. The particularity of this method is to use casein as the protein 


kp 
Pd y Q) So PEN 
Q FF / \ emulsifying 
XB ¢ + ( ) 
= | S. < 
Casein Oil droplet of 
molecules monomer/alkane 


washing 


Clean 
microcapsule 


multilayer to protect the microparticles from re-coalescencing 
before MPCM is finally stabilized. The fabrication process is shown 
in Fig. 2. They used n-hexadecane as the core material. It was 
shown that the casein was very effective in protecting the 
agglomeration of emulsion. Further measurements indicated that 
the melting temperature of MPCM was increased and latent heat of 
MPCM was smaller than the bulk material which they attributed to 
the presence of the impurity. The decrease of the latent heat was 
always a drawback of MPCM because of the presence of the shell. 

Zou et al. [8] fabricated the MPCM by the interfacial 
polycondensation using hexadecane as the core material. It was 
found that the average diameter of MPCM was about 2.5 wm, and 
the latent heat was drastically decreased to 66.09 J/g compared to 
236 J/g of pure hexadecane; while the melting temperature was 
not seriously affected. 

MPCMs are generally synthesized with melamine-formalde- 
hyde (MF), urea—formaldehyde (UF) or gelatin—formaldehyde as 
the shell materials. However, formaldehyde is harmful to the 
health, the remnant of which should be reduced to meet the 
limitation by the standard when MPCM is used for textile. Li et al. 
[9] prepared the melamine-formaldehyde (MF) prepolymer by 
incorporating formaldehyde once and melamine three times. Urea 
was used to remove the remnant of formaldehyde after the in situ 
polymerization finished. The morphology of the fabricated MPCM 
was investigated for different dropping rate of MF prepolymer 
into the n-octadecane emulsion. The characterization showed 
that the average diameter was about 2.2 wm with very narrow 
range and the MPCM contained 59 wt% of n-octadecane with the 
enthalpy of 144 J/g. The measurement showed that the remnant of 
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Fig. 2. A simple method for making MPCM proposed by Ai et al. [7]. 
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formaldehyde was only 67.1 mg/kg, which implied that MPCM 
had the potential application for the textile. 

In addition to the paraffin used as the core material, coco fatty 
acid mixture can also be the candidate for MPCM. Ozonur et al. [10] 
used inexpensive natural coco fatty with the melting temperature 
in the range of 22-24 °C as the core material and several kinds of 
possible shell materials were tested by the coacervation techni- 
ques. Results indicated that gelatin-gum Arabic was the best shell 
material in microencapsulation, and the diameters of the 
fabricated MPCMs were around 1000 um. The MPCM would be 
very cheap by using the proposed cheap core and shell materials, 
which opened a new way to reduce the cost of MPCM. 

A lot of effort has been devoted to the determination of the 
thermal properties of MPCM which are obviously very much 
different from those of the bulk materials. In the energy storage 
and transportation applications, the thermal properties including 
latent heat, subcooling degree, melting temperature, viscosity, etc. 
are the most important ones. 

Subcooling is always a serious problem in PCM research and 
application. One of the drawbacks of MPCM is the possibility of 
increasing the subcooling degree, which should be effectively 
suppressed before MPCM can be used as the energy storage and 
transportation medium, because the subcooling will drastically 
deteriorate the system performance and reduce the energy 
efficiency, for example, increasing the refrigerating load and 
reducing the efficiency of the refrigerating system for chilling 
MPCM in air conditioning applications. 

Alvarado et al. [11] experimentally characterized the thermal 
behavior of the bulk and microencapsulated n-tetradecane by 
using DSC. The subcooling degree was specially studied and 
compared, and several methods were tried to reduce the 
subcooling degree. Silica fume as the nucleating agent was test 
and the results showed that its presence could significantly 
suppress the subcooling for the bulk material, while 0.2 wt% of 
silica fume was not effective for MPCM, which indicated that a 
suitable nucleating agent is necessary for the future study. A 
comparative study by using a more stable nucleating agent- 
tetradecanol with 1-4 wt% was carried out, and it was showed that 


Table 3 
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2 wt% of tetradecanol was already very effective in suppressing 
subcooling. However, the mass fraction of the nucleating agent 
should not be too high because it will reduce the latent heat 
significantly with increasing the concentration. 

There are several kinds of materials which can be used as the 
nucleating agents to suppress the subcooling of MPCM. Fan et al. 
[12] prepared three kinds of MPCMs (n-octadecane as the core 
material) with sodium chloride, 1-octadecanol and paraffin as the 
nucleating agents, respectively. It was revealed that the subcooling 
could be effectively suppressed by increasing the concentration of 
the nucleating agent to 6, 9 and 20 wt% for sodium chloride, 1- 
octadecanol and paraffin, respectively. However, addition of 
sodium chloride and 1-octadecanol as the nucleating agents had 
the drawbacks on the morphology and dispersibility of MPCM, 
which resulted in rough shell surface and agglomeration. Never- 
theless, addition of paraffin to the core material had little effect on 
the morphology and dispersibility of MPCM. This point may be 
very important because rough surface of MPCM will increase the 
circulating pumping power which is undesirable to the practical 
application for the energy storage and transportation. 

Table 3 shows the summary of the core and shell materials, 
nucleating agents and the size distribution of MPCM. Some 
information is from the commercial available products, and the 
most are the laboratory-scale test samples. 


2.2. Fundamental features of MPCS 


2.2.1. Thermal properties of MPCS 

When the MPCM is dispersed into the carrier fluid, e.g. water, a 
kind of suspension named as microencapsulated phase change 
slurry (MPCS) is formed. Due to the relatively large surface area to 
volume of MPCM, better heat transfer performance can be 
achieved. Thus, MPCS can serve as both the energy storage and 
heat transfer media, in addition, the agglomeration of MPCM can 
be avoided due to the encapsulation. It is obvious that, however, 
the heat transfer and fluid flow characteristics are very important 
for the MPCS system design. As a kind of suspension, the thermal 
properties of MPCS are different from those of the bulk PCM 


Summary of the core and shell materials, nucleating agents and size distribution of MPCM. 


Reference Core material Shell material 


Nucleating agent Size distribution (wm) Comments 


Alkan et al. [5] Docosane Polymethylmethacrylate 


(PMMA) 


Jin et al. [6] Paraffin wax Urea-formaldehyde 

Ai et al. [7] n-Hexadecane Polystyrene 

Zou et al. [8] Hexadecane Shell was formed by 
interfacial reaction 
between TDI and EDA 

Li et al. [9] n-Octadecane Melamine-formaldehyde 
(MF) 

Ozonur et al. [10] Coco fatty Gelatin-gum Arabic 

Alvarado et al. [11] n-Tetradecane gelatin 


Fan et al. [12] n-Octadecane 


Zhang et al. [13] Urea-melamine- 
formaldehyde 


Amino plastics 


8.8 vol% cyclohexane 

and 91.2 vol% n-octadecane 

Wang et al. [16] 1-Bromohexadecane (C;¢H33Br) 

Inaba et al. [20] 

BASF (Micronal™ 
DS 5000, 5007) 

Microtek 


n-Docasane, n-tetradecane 
Paraffin wax PMMA 


Paraffin wax Polymer 


0.16 (average) 


5-20 (depending on 
synthesis condition) 


3-15 Latent heat was drastically 
decreased due to the impurity 
1.26-10 Latent heat was drastically 


(2.5 average) decreased 


2.2 (average) Formaldehyde was reduced 


1000 
90-125, 
(100 average) 


Very cheap 

0.2 wt% of silica fume 
(not effective for MPCM) 
tetradecanol with 1-4 wt% 
(effective) 


0.2 wt% silica fume, 
2 wt% tetradecanol, etc. 


Sodium chloride, 
1-octadecanol and 
paraffin 
0.4-5.6 


2.0-31.7 

(10.1 average) 

15, 1.5 

5 Commercially available 


Commercially available 
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materials and carrier fluids, which are very essential for the 
evaluation of the flow and heat transfer characteristics of MPCS. In 
the following context, brief discussion of the determination of the 
thermal properties of MPCS is first introduced, which mainly 
include the thermal conductivity, viscosity and specific heat. 

In general, the thermal conductivity of the static MPCS with low 
volume fraction can be evaluated by Maxwell’s equation for- 
mulated as 


2kp + kp + 2h(kp — ke) 
2k + kp _ (kp = kr) 


ks = kç x (1) 
where k, and kç are the thermal conductivities of the bulk slurry 
and carrier fluid, and @ is the volume fraction of MPCM. The 
thermal conductivity of the MPCM particle kp can be estimated by a 
correlation describing composite sphere feature [14]. Under the 
flowing condition, the effective thermal conductivity is generally 
higher than that predicted by Maxwell’s equation due to the 
interaction between the particle and fluid. Then, a general 
correlation summarized in reference [15] can be applied depend- 
ing on the particle Peclet number 

ke 
— = 1 + BoPep (2) 
ks 

where ke is the effective thermal conductivity of the PCS, Pep = 
é. d/a is the particle Peclet number, ë is the shear rate, a; is the 
thermal diffusivity, and the values for B and m depend on the 
particle Peclet number. In general, the thermal conductivity of 
MPCM particles is not very large and needs to be enhanced. It has 
been proposed that a certain amount of nano-particles are added to 
PCM so as to enhance the thermal conductivity of MPCM as in the 
case of nano-fluids. 

The viscosity of MPCS is critical in determining the pressure 
drop when MPCS is used to transport energy. When the volume 
fraction of MPCM is less than 25%, the MPCS is generally considered 
as the Newtonian fluid. The viscosity is generally measured by 
using viscometer, and the results [16,17] indicated that the MPCS 
showed the Newtonian fluid feature and the apparent viscosity 
was about 1.5-10 times of the water as the mass fraction of MPCM 
increased from about 5 to 30%. The MPCS can be treated as the 
homogeneous fluid at low volume fraction, and the apparent 
viscosity of the slurry can be estimated by the following equation 
[18] which includes the effect of the interactions between MPCM 
particles, carrier fluid and even the tube wall. 

Hs (1 @-ag?) (3) 
Me 

where A is a parameter depending on the shape and rigidity of the 
particles. The pressure drop can be calculated by the equation for 
the single phase flow by using the apparent viscosity, however, the 
effect of non-Newtonian rheology will be more evident as the 
volume fraction increases beyond 30%. It has been demonstrated 
that MPCS shows Newtonian rheology in the most studied volume 
fraction range, nevertheless, it displays mostly non-Newtonian 
rheology for clathrate hydrate slurry even at low volume fraction, 
which will be discussed in detail in Section 3. 

The specific heat of the MPCS without phase change can be 
simply estimated from the mass and energy balance approach 
which can be formulated as 


Cp s = WpcemCpem + Wshell she + WC (4) 


where Wpcm, Wsheil and wp are the mass fractions of the PCM, shell 
and carrier fluid, respectively, and Wpcm + Wspe + Wp = 1. How- 
ever, when the phase change is involved in the estimation of the 
heat transfer performance, the effect of the latent heat should be 
properly treated and an effective specific heat is used which can be 


written as if all the MPCM has been melted 
Ce = Cop + Qa 
where AT is the phase change temperature range. 


2.2.2. Flow and heat transfer characteristics of MPCS 

Charunyakorn et al. [15] proposed a theoretical model to study 
the forced convection heat transfer of MPCS, in which several 
dimensionless parameters, like the bulk Stephan number, particle 
volume fraction, modified Peclet number and so on, are included. 
The results showed that the bulk Stephan number and particle 
volume fraction were the most important parameters which 
affected the heat transfer of MPCS. 

Goel et al. [19] conducted experiments to investigate the 
laminar forced convection heat transfer of MPCS in order to 
provide the useful information which could be directly used by 
those who were not involved in this field. The Reynolds numbers 
studied were about 200 and 1000. The experimental results 
indicated that about 50% reduction in wall temperature could be 
realized by using MPCS compared with single phase fluid flow and 
the Stephan number, defined in the following equation, was the 
most dominant parameter in the convection heat transfer similar 
to the results of Charunyakorn et al. [15]. 


= CpslqwR/ks| 
pL Pp/Ps 


Due to the temperature control application purpose, either heat 
transfer coefficient or Nusselt Number was not chosen to 
characterize the heat transfer performance. A dimensionless wall 
temperature defined as Owx=(Twx—Tji)ks/(qwD) was used to 
represent the wall temperature. Shown in Fig. 3 is one typical 
result under different Stephan numbers when the volume fraction 
of MPCM is 10% and diameter is about 100 um. Although the 
comparison of the experimental results to theoretical ones of 
Charunyakorn et al. [15] was made, it was found that the 
agreement was only qualitatively well due to the complexity of 
the problem. The effects of the other influential parameters, like 
volume fraction, diameter of MPCM and the homogeneity of MPCM 
particles, have also been investigated and discussed. 

Inaba et al. [20] studied the heat transfer characteristics of 
MPCS composed of different sized particles in a circular tube. 
MPCM with a nominal size of 15 wm (n-docasane as the core 
material) in diameter was blended with smaller sized MPCM of 
1.5 wm (n-tetradecane as the core material) in diameter. The flow 
of MPCS exhibited smaller friction factor than that of the pure 
water as the drag-reduction additives which could reduce the flow 
resistance was added. The MPCS with different sized particles 
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Fig. 3. Comparison of the experimental and theoretical results for the dimensionless 
wall temperature [19], volume fraction 10%, Re = 200, d= 100 pm. 
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showed better heat transfer performance than that with mono- 
sized particle, which was due to the microconvection around the 
particles. 

The advantage of application of MPCS to the energy storage and 
transportation is that the phase change with latent heat is involved 
and thus the effective specific heat of the fluid is remarkably 
increased which results in the heat transfer enhancement. 
However, one has to take the pumping power into consideration 
when putting MPCS into the practical application. Inaba et al. [20] 
found that the flow drag reduction could be achieved in the 
turbulent flow regime, while the pressure drop was increased in 
the laminar flow regime in the case of higher mass fraction of larger 
sized particles. Thus, the ratio between the heat Q transported by 
MPCS (including both latent heat and sensible heat) and the 
pumping power W of the circulating pump has to be investigated. 
They showed experimentally that Q/W for MPCS was larger than 
that for water in the turbulent flow regime, while in the laminar 
flow regime, Q/W decreased from larger than that for water to less 
than that for water as increasing the fraction of the larger MPCM 
particles, as shown in Fig. 4. This fact indicated that both the energy 
transportation and pumping power should be simultaneously 
considered in the application of MPCS. 

MPCM should be endurable before it can be put into long-term 
operation, because it might be broken by the circulating pump. 
Alvarado [21] showed that MPCM (n-tetradecane as the core 
material) with the diameter of smaller than 10 um could be 
endurable and impact-resistant by using cavity pump to do the 
durability test. The pressure drop and heat transfer characteristics 
were also experimentally measured. It was showed that the 
pressure drop of MPCS at low mass fraction was much smaller than 
that of water, while the differences decreased as the mass fraction 
increased, which might imply that there was a possible drag 
reduction effect of MPCS. The heat transfer coefficient of MPCS 
reached its peak when the temperature was close to the melting 
temperature of n-tetradecane, which indicated that additional 
specific heat caused by the phase change enhanced the heat 
transfer. However, the heat transfer coefficient of MPCS was 
generally smaller than that of water under the same condition, 
which might be attributed to the turbulence reduction and smaller 
thermal conductivity of MPCS. The using of the enhanced tube 
might be effective in enhancing the heat transfer. 

Wang et al. [16] presented a comprehensive study of the heat 
transfer of MPCS in circular tube. They used 1-bromohexadecane 
(Cie6H33Br) as the core material and amino plastics as the shell 
material for microencapsulated particles. The experimental heat 
transfer coefficient was compared with Gnielinski correlation for 
single phase flow and it was found that the correlation for single 
phase flow could not properly describe the heat transfer of MPCS. 
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Fig. 4. Variation of the ratio of the heat Q transported by MPCS to the pumping 
power W with Re’ [20], where Re’ is the modified Reynolds number based on power 
viscosity law, C in the figure is the mass fraction. 


Although there were several empirical heat transfer correlations 
for solid-liquid two-phase slurry flow reported in the literature, 
however, they were only applicable to the heat transfer of slurry 
without phase change. In the case of MPCS, the authors finally 
proposed their own correlations based on the experimental data by 
including the Stephan number. The comparison of the experi- 
mental data and correlation showed that the relative errors for the 
laminar and turbulent heat transfer were within +10%. It was 
concluded that the heat transfer was enhanced for MPCS and the 
average Nusselt number was about 1-2.5 times higher than for the 
pure water in the case of the turbulent flow. 

The phase change problem of MPCM attracted quite a lot of 
attention because of its fundamental importance. Zhao et al. [22] 
carried out a parametric analysis of the heat transfer enhancement 
of the laminar flow of MPCS in the circular tube under the constant 
wall temperature condition. They adopted the definition of the 
modified local heat transfer coefficient and modified Nusselt 
number to avoid unreasonable inflexion of the local heat transfer 
coefficient and Nusselt number. Totally, six factors including the 
Stefan number, volume fraction of MPCM particles, ratio of the 
particle radius to tube radius, subcooling, dimensionless phase 
change temperature range and the bulk flow Reynolds number 
were taken into consideration in the analysis. The analytical results 
indicated that the Stephan number, volume fraction and size of the 
MPCM particle were the most important influential factors for the 
heat transfer enhancement, and the enhancement up to 4 times 
could be achieved by employing MPCS as the working fluid. 

MPCM can be used to increase the energy storage density of the 
heat storage system. A mathematical model was developed to 
study the thermal behavior of a MPCS latent heat storage system 
on the basis of the energy conservation and heat transfer [23]. It 
was revealed that the convective heat transfer coefficient between 
the MPCM particles and the carrier fluid (water) played a very 
critical role in the charging and discharging cycles. Some other 
influential factors including MPCM mass fraction and particle 
diameter were also investigated. These results were very helpful in 
optimizing the thermal system although some simplifications 
were made in the analysis. 

Yamagishi et al. [24] carried out quite detailed tests of the 
MPCM particles with the core material of n-tetradecane and n- 
dodecane (with the melting temperatures at 5.0 and —13.5 °C as 
reported by the authors, respectively). Although the above two 
bulk materials did not have evident subcooling degrees, it was 
found from the test results that the subcooling was quite 
significant as the diameter of the particle decreased, which could 
be even as large as up to 15°C, as shown in Fig. 5. The 
heterogeneous nucleation is generally preferred than the homo- 
geneous nucleation because the former occurs at higher tempera- 
ture. In order to solve the subcooling problem, the nucleating agent 
of 1-tetradecanol for n-tetradecane was applied, which effectively 
decreased the subcooling and reduced the cooling time. The 
tolerance against the shear force induced by pumping is very 
essential for practical application. Breakage tests were carried out 
for different sized particles, which revealed that the smaller 
particles could survive the test. MPCS with high volume fraction 
generally has quite large apparent viscosity which consequently 
increases the pressure drop and pumping power. However, the 
drag-reduction surfactants can be used to reduce the pressure loss 
in the tube. Shown in Fig. 6 is the variation of the friction factor 
with the Reynolds number for MPCS with the surfactant, which 
evidently shows that the flow characteristics display the similar 
feature to that of the Newtonian fluid. 

It can be seen that it is very useful to disperse MPCM into the 
carrier fluids to form MPCS for the thermal storage and manage- 
ment application. The convective heat transfer coefficient is 
generally larger than that of the single phase fluid flow because 
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Fig. 5. Variation of the melting and crystallization temperatures with the size of 
MPCM [24]. 
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Fig. 6. Variation of the Darcy friction factor with the Reynolds number [24]. 


of the increase of the specific heat by the phase change. The 
pressure drop can be generally predicted by the single phase flow 
correlation when the volume fraction of MPCM is not very large, 
and the pressure drop can be effectively reduced by using the drag- 
reduction additives in the case of high volume fraction of MPCM. 


2.3. Application of MPCS to the thermal systems 


The application of MPCS to the thermal systems mainly 
includes its using in the thermal storage and transportation 
system, and in the thermal controlling aspects and so on. One 
typical example is the secondary loop for refrigeration and air 
conditioning. 

An assessment of the performance of the embedding MPCS 
energy storage and evaporative cooling into a cooled ceiling 
system was carried out by Wang et al. [25]. The core material of 
MPCM was hexadecane with a melting temperature of 18.1 °C, 
which enabled the possible integration into the cooled ceiling 
system usually working at 14-20 °C. In the research, MPCS was 
only used for the energy storage, but not used as the working fluid 
for the secondary loop. The results indicated that the annual 


energy saving potential up to 80% could be achieved in north- 
western China. 

The comparative studies of the energy charging and releasing 
performances of PCM in an annular space were carried out by 
Balikowski and Mollendort [26] by using Climsel 28 (ordinary 
PCM) and microencapsulated Thermasorb 83. It was found that the 
heat transfer rate of latter was about 25% higher in the charging 
period and also 20% higher in the discharging period, which 
indicated that MPCM could significantly increase the heat storage/ 
retrieval rates. 

Griffiths and Eames [27] have conducted experimental tests by 
using MPCS (products of BASF) as both the energy storage and 
transportation media in a cooled ceiling system. It was concluded 
from four months-long experiments that 40 wt% MPCS could 
effectively reduce the flow rate while maintaining the constant 
cooling temperature, thus the total pumping power was reduced 
compared to the large flow volume of the water. 

Fossett et al. [28] has proposed to use the MPCM embedded 
heat sink instead of a solid aluminum plate as the passive heat sink 
for the avionic application. The advantage of using MPCM in the 
heat sink over the bulk PCM lies in eliminating the volume change 
during melting or solidification. Both the experimental and 
analytical investigations on a heat sink with a size of 
12.7 cm x 22.9cm in basement and 2.54cm in thickness were 
carried out, and 270 g MPCM with the melting temperature close 
to the temperature range of the cooling application was used. The 
heating load of 90 W was applied to the heat sink and was lasting 
for 10 min. It was found that the temperature of the MPCM heat 
sink could be well maintained below the safety temperature of the 
electronic device at 71 °C in 10 min at an initial temperature of 
52 °C which was about 10 °C smaller than the melting temperature 
of MPCM. However, it was only one minute when the temperature 
of the heat sink without MPCM reached 71°C under the same 
initial condition. The analytical results agreed quite well with the 
experimental ones. This technology is very useful for the design of 
the light airborne heat sink. 

A quite large-scale application of MPCS to cooling in Narita 
airport of Tokyo was described by Shibutani [29], which was 
implemented by the co-operation between Narita airport and 
Mitsubishi Heavy Industries Ltd. The initiation of the project was 
motivated by the phase-out of refrigerants R11 and R22. If the 
system was re-charged with R134a and R123, there would be the 
decrease of the cooling capacity and performance without 
replacing the chiller system. Thus, there was insufficient cooling 
power supply in the daytime. However, this problem could be 
overcome by incorporating the cold storage system into the 
existing chiller system. The temperatures of the chilled water in 
the supply and return pipes of the existing system were 5 and 
12 °C, respectively. The MPCM with the average diameter of 2 ym 
and phase change temperature between 7 and 8 °C was therefore 
selected for the cold storage system. Shown in Fig. 7 is a typical 
temperature variation of the MPCM used for the cold storage and 
cold release, in which an evident phase change part could be easily 
observed and the significant increase of the cold storage density 
was achieved compared to the chilled water. A storage tank with 
the volume of 970 m? was installed onsite and was coupled to the 
chiller through a plate heat exchanger. The system operation 
showed that the cold storage density was 67 MJ/m? with MPCS, 
larger than 21 MJ/m? with cold water, although much smaller than 
167 MJ/m? with ice. While the energy efficiency and running cost 
for the former two cases were similar, better than the ice storage. 

There are already commercial available MPCMs in the market 
with the typical products from BASF and Microtek company. 
Gschwander et al. [30] dispersed BASF MPCMs, as shown in Table 3, 
into water to form MPCS and studied the durability after pumping 
for several weeks at a rate of 800 cycles/day. The long-term tests 
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Fig. 7. Temperature variation of the MPCM used for the cold storage and cold release 
[29]. 


indicated that the centrifugal pump was suitable for pumping 
MPCS and caused less damage to the shell of MPCS. Their initial test 
indicated that MPCM might be broken when transported through 
pump, while this could be tackled by increasing the thickness of 
the shell and reducing the diameter, although it resulted in less 
melting enthalpy of MPCS. The pressure drop and temperatures 
were recorded in the test in order to evaluate the performance of 
MPCS, and the thermal tests implied that MPCS was especially 
advantageous for the applications with a small temperature 
difference. 

Many investigations have demonstrated that MPCM can be 
readily dispersed into carrier fluids to form MPCS, which can be 
used as both the heat transfer and energy storage media. 
Nevertheless, much effort is still needed to improve the 
performance, such as specific heat, thermal conductivity and 
durability, and to understand the flow and heat transfer 
characteristics of MPCS, so as to promote the applicability at the 
industrial scale. 


3. Clathrate hydrate slurry (CHS) 


The other big category of the PCS is clathrate hydrate slurry 
(CHS). Although clathrate hydrate itself has already been subjected 
to the intensive study in the field of chemistry, its application to 
the energy storage and transportation is still very new and at its 
early stage. Clathrate hydrate is generally produced under high 
pressure and low temperature condition, however, some semi- 
clathrate hydrate, for example, semi-clathrate hydrates of 
quaternary ammonium salts (e.g. tetrabutylammonium bromide, 
TBAB in abbreviation) can be formed under the condition of one 
atmosphere pressure and room temperature. Due to the fact that a 
phase change with the latent heat up to several hundred kJ/kg is 
associated with the formation and dissociation of the semi- 
clathrate hydrate, it is possible to be employed in the energy 
storage and transportation application. For example, the cold 
storage capacity of TBAB CHS was about 2-4 times larger than that 
of the cold water in the temperature range of 5-12 °C depending on 
the mass fraction of the hydrate crystals in the slurry, which 
implies that the flow rate can be decreased and therefore the 
pumping power can also be reduced. Shown in Fig. 8 is the 
comparison of the enthalpy changes of the TBAB CHS, TBAC CHS 
and water in the temperature range of the air conditioning 
application, i.e., 5-12°C. The enthalpy changes of CHS are 
calculated from the DSC curves and it is obvious that the enthalpy 
change of CHS is much larger than that of the water. 

Although the research of the semi-clathrate hydrates of 
quaternary ammonium salts can be traced back to several decades 
ago, it does not draw the intensive attention until recent years. This 
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Fig. 8. Comparison of the enthalpy changes of the TBAB CHS, TBAC CHS and water in 
the temperature range of 5-12 °C, the base point is taken at 12 °C. 


section provides the review of the studies of CHS used as PCS, 
focusing on the fundamental properties of CHS, the generation, 
flow and heat transfer characteristics and the applications. 


3.1. Fundamental properties of clathrate hydrate crystals 


3.1.1. Structures and shapes of hydrate crystals 

Clathrate hydrates of tetraalkylammonium salts with simple 
anions (such as halides, sulfate, formate, etc.) were firstly reported 
by Fowler et al. in 1940 [31]. Later, many researchers carried out a lot 
of investigations on clathrate hydrates, most of which were focused 
on the crystal data and structure refinement from chemical aspects. 
Nonetheless, some thermal properties were also reported. 

Clathrate hydrate is a chemical compound in which the water 
molecules (“host”) form lattice structures and molecules of another 
substance (“guest”) fill the lattices. It is known that salt anions are 
bound to water molecules through hydrogen bonds and form a 
water-anion framework in the semi-clathrate hydrate [32]. In the 
following context, TBAB is mainly discussed as an example because 
its clathrate hydrate was subjected to the most intensive studies. 
Even under one atmosphere pressure and room temperature 
condition, TBAB in water can form a semi-clathrate hydrate in 
which TBAB molecules are both the host and guest, as shown in Fig. 9. 
There are still small dodecahedral empty cages in TBAB clathrate 
hydrate. Hence, TBAB clathrate hydrate is able to uptake, store and 
separate small gas molecules under different thermal conditions, 
which have spawned a lot of researches on this topic recently [33,34]. 


Fig. 9. Molecular structure of TBAB semi-clathrate hydrate. 
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Fig. 10. Two types of TBAB hydrate crystals [35]. 
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Fig. 11. Two temperature ascents in the cooling down of TBAB aqueous solution. 


TBAB CHS appears white and is a kind of solid-liquid 
suspension like ice slurry. There are generally two different types 
of TBAB hydrate crystals, as shown in Fig. 10, and their existences 
have been confirmed by the facts that the temperature ascends 
twice during the cooling down of TBAB aqueous solution. Shown in 
Fig. 11 is the cooling curve of TBAB aqueous solution with two 
temperature spikes which indicate two different phase changes 
corresponding to the formation of two different hydrate crystals 
with different transmittances, refractive indices and crystal 
morphologies. Type A is columnar shape and type B is an irregular 


Type A TBAF hydrate crystals 


form composed of thin crystals [35]. Shimada et al. [36] observed 
the growth of TBAB hydrate crystals using a digital high definition 
microscope, and showed that the growth of TBAB hydrate was 
highly dependent on the interfacial kinetics. More information 
about structures and shapes of TBAB clathrate hydrate crystals, etc. 
can be found in the references [32,37]. Similar to TBAB, Rodionova 
and Manakov [38] also observed two types of clathrate hydrates of 
tetrabutulammonium fluoride (TBAF), and they isolated both 
hydrate crystals from the aqueous solution and their appearances 
are shown in Fig. 12 with tetragonal structure and cubic structure, 
respectively. 


3.1.2. Phase diagram of CHS 

Through cooling down the salt aqueous solution with strong 
agitation, hydrate crystals will nucleate in the liquid. Hydrate 
crystals do not appear until the temperature reaches a certain 
degrees lower than the crystallization temperature due to the 
subcooling phenomenon. Once the nucleation occurs in the 
subcooled aqueous solution, the subcooled state will gradually 
disappear due to the release of the latent heat in the exothermic 
process, and the mixture temperature will ascend to the phase 
change temperature of the original aqueous solution. If the original 
salt concentration of the solution is equal to the salt mass fraction 
in hydrate crystal (e.g. 40.8 wt% for TBAB-26H20), the salt 
concentration in liquid phase will remain constant even after 
the hydrate crystals form. Then the hydrate crystals will keep 
forming until all solution is crystallized, therefore the CHS 
temperature remains nearly constant due to the latent heat 
released. Based on this phenomenon, the congruent concentration, 
the corresponding melting temperature and the hydration number 
can be determined. Other melting temperatures at incongruent 
concentrations can also be determined from the cooling curve. 
However, after the hydrate crystals form, the salt concentration in 
liquid phase changes and lower temperature is needed to sustain 
the crystallization, so the CHS temperature can not be kept 
constant but decreases due to the cooling. 

The heating in DSC test was used by many researchers to 
determine the melting temperatures and latent heats of hydrate 
crystals while the cooling was not suitable due to the subcooling 
phenomenon. Typical DCS curves are shown in Fig. 13. One of the 
important factors that affect the measurement accuracy is the 
heating rate which has been investigated in many studies. 
Nakayama [37] and Rodionova and Manakov [38] used the heating 
rate of 0.5 °C/min in the measurement, while Yamazaki et al. [39] 
applied the heating rates of 2 °C/min and 0.1 °C/min in the case of 
low salt concentration. It is generally accepted that higher heating 
rate results in a better measuring sensitivity, but the penalty is that 
the resolution will be decreased. 

According to Oyama et al. [35], DSC was not suitable for 
determining the phase diagram of TBAB clathrate hydrates, 
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Fig. 12. Two types of TBAF hydrates crystals [38]. 
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Fig. 13. Typical DSC heating curves of CHS [37], X in the figure is molar fraction. 


because there are two kinds of TBAB hydrate crystals and their 
melting temperature peaks are too close to be separated. They 
measured the melting temperatures of these two types of hydrate 
crystals by two different techniques. For type A hydrate, the 
common method was used and a heating rate of 2 °C/h was applied 
for CHS; for type B hydrate, they grew and separated the crystals 
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and raised the temperature at a step of 0.1 °C, and observed the 
crystals by a digital CCD microscope to determine the dissolving 
temperature, which was indicated by the surface of crystals 
becoming visibly rough. 

Fig. 14 shows the phase diagrams of four typical clathrate 
hydrates which have been studied as the thermal energy storage 
media except tetrabutylammonium chloride (TBAC). A lot of 
information about CHS can be obtained from the phase diagram. 
Salt mass fraction in the liquid phase x; can be determined from the 
corresponding CHS temperature T on the phase diagram, then 
hydrate mass fraction x, and volume fraction @ can be calculated 
by the following equations: 


Xo — XI 


XH 
XN — XI 


(7) 


_ XH/ Pu 
Xu/Py + (1 — Xn) /A 


(8) 


where xo is the initial salt mass concentration in the aqueous 
solution, xn is the salt mass fraction in the hydrate crystal, p4 and pı 
are the densities of the hydrate crystal and liquid, respectively. 


3.1.3. Thermal properties of hydrate crystals and CHS 

Thermal properties of the hydrate crystal are essential for 
understanding the fundamental features of clathrate hydrate and 
further application. Given in Table 4 are the thermal properties of 
TBAB, TBAF, TBAC and TME hydrate crystals, including the 
hydration number, melting temperature, density, latent heat, 
specific heat and thermal conductivity. 

The latent heat and specific heat of the clathrate hydrate are 
usually determined from the DSC curves, in which the effect of the 
ice should be carefully considered. Nakayama [37] analyzed the 
cooling curves of TBAC aqueous solution measured by DSC, and the 
hold point which was after the formation of hydrate and before the 
formation of ice was determined to separate the enthalpy change 
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Fig. 14. Phase diagrams of four CHS (a) TBAB hydrate, hydrate number Ntype_q = 26, Ntype_p = 36 (38 according to Oyama [35]); (b) TBAC hydrate, N = 30 [37]; (c) TBAF hydrate, 
N = 30 according to Nakayama [37], N = 29 according to Dyadin [40], Ntype_a = 30 and Niype_s = 33 according to Rodionova [38]; (d) TME hydrate, N = 3. 
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Table 4 
Characteristics of TBAB/TBAC/TBAF/TME hydrate crystals. 
Authors Hydration Melting Density Latent heat Specific heat Thermal 
number temperature (°C) (kg/m?) (kJ/kg) (kJ/(kg K)) conductivity 
(W/(mK)) 
TBAB Lipkowski et al. [32] 24/26/32/36/2/3 12.4/12.2/11.7/9.5/-/- 
Oyama et al. [35] 26/38 12.0/9.9 193.2/199.6 *1.86-2.61/2.0-2.54 
Darbouret et al. [48] 26/36 12.3/9.6 1082/1067 
Hayashi et al. [49] 
Ogoshi and Takao [55] 26/36 11.8/- 1082/1030 193/205 2.22/- 0.42/- 
TBAC Nakayama [37] 30 15.0 200.7 
Dyadin and Udachin [40] 29.7/32.1 15.0/14.7 1034/1029 
TBAF Nakayama [37] 30 28.3 229.7 
Rodionova and Manakov [38] 29.7/32.8 222 223.1/240.5 
Dyadin and Udachin [40] 28.6/32.3 27.4/27.2 1057/1035 
TME Yamazaki et al. [39] 3 30.0 1120/1090 218 °2.75/3.58 0.65/0.21 


* Between —20 and 0.2°C. 
> At 10 and 50°C. 
© At 22 and 56°C, respectively. 


for ice and hydrate, as shown in Fig. 15. The author obtained the 
total enthalpy change AH, and the hydrate enthalpy change AH, 
then the latent heat L and hydration number N can be calculated 
by: 


E AH, 
L= ia (AH; — MA) [Lice (3) 
N= (Mpo — (AH; — AH2)/Lice)/My,0 (10) 


Msait /Msalt 


Another way of using DSC to measure the latent heat of the 
hydrate crystal directly was proposed by Oyama et al. [35]. The 
authors separated the two types of TBAB hydrate crystals from the 
aqueous solution, and these crystals were further stabilized by 
quenching in liquid nitrogen and then were measured by using 
DSC directly. Therefore the sole latent heat without the inter- 
ference of the enthalpy change of solution was obtained. The key 
point of this method is the purity of the hydrate crystals because 
ice would more or less adhere to the hydrate crystals. Hydrate 
crystals growing at their congruent concentration will be in high 
purity, however, not all hydrates are congruently crystallized, such 
as type B TBAB hydrate, etc. Oyama et al. took the influence of the 
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Fig. 15. The method used to determine the latent heat of hydrate on DSC curve [37]. 


ice remnant on the surface of type B TBAB hydrate crystals into 
consideration, and obtained the final results of the latent heats for 
TBAB hydrate crystals. Moreover, the specific heats of TBAB 
hydrate crystals were determined from the same DSC results. It 
was considered that the larger latent heat and specific heat of type 
B TBAB clathrate hydrate are the results of larger hydration 
number. 

Thermal conductivity of CHS was rarely reported in the 
literature, however, it can be calculated by Maxwell’s equation 
shown in Eq. (1). And the approaches depicted in the sections for 
determining the thermal conductivity of MPCS can also be 
applicable to CHS. 


3.1.4. TBAB clathrate hydrate combined with gas molecules 

Semi-clathrate hydrate is able to uptake, store and separate low 
molecular weight gases, which opens a new way for the gas storage 
and transport. Arjmandi et al. [41] studied the equilibrium data of 
hydrogen, methane, nitrogen, carbon dioxide and natural gas in 
TBAB clathrate hydrate, which showed that these double hydrates 
were more stable than the single TBAB hydrate. The dissociation 
enthalpy of gas hydrate is generally larger than that of salt 
clathrate hydrate, for example, 374 kJ/kg for CO, hydrate. Usually, 
the equilibrium pressure of gas hydrate is very high. However, 
tetrahydrofuran (THF) or TBAB can be added to reduce the 
equilibrium pressure of gas hydrate, as shown in Fig. 16, which is 
obviously more favorable for application. 

Delahaye et al. [42] investigated the effects of THF on the 
equilibrium pressure decrease and dissociation enthalpy enhance- 
ment of CO, hydrate. And 3.8-15 wt% THF was added to water, 
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Fig. 16. Equilibrium pressure reduction of CO2 hydrate with TBAB/THF in water [42- 
44]. 
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which showed that the equilibrium pressure of CO2 hydrate has 
been reduced by about 80% and the dissociation enthalpy of the 
THF + CO2 double hydrate in the unit of kJ/mol was nearly two 
times of that of single CO2 hydrate. 

Lin et al. [43] applied the differential thermal analysis (DTA) to 
measure the dissociation enthalpy of TBAB + CO2 double hydrates. 
The composition of hydrate formed with 9.01wt% TBAB + CO, is 
2.51CO2:TBAB-38H20 and its dissociation enthalpy is 313.2 kJ/kg, 
smaller than that of the single CO, hydrate, but larger than that of 
TBAB clathrate hydrate (about 200 kJ/kg). Li et al. [44] tested the 
equilibrium data of TBAB solution combined with CH4, which 
showed that the equilibrium pressure was reduced drastically by 
52-96% compared to the single CH, hydrate and the equilibrium 
temperature was increased by 4-17 °C compared to that of TBAB 
clathrate hydrate. 


3.2. Generation of clathrate hydrate crystals 


The basic steps to generate the hydrate crystals from the 
aqueous solution include supersaturation, crystals nucleation 
and growth. Following these three steps, attrition, agglomeration 
and ripening are needed to make the hydrate particles to take the 
prescribed shape and size. Hydrate can be formed only if there is 
a sufficient driving force, therefore supersaturation state is 
necessary. Usually, supersaturation can be achieved by decreas- 
ing the temperature to lower than the crystallization tempera- 
ture or by changing the pressure to shift the equilibrium 
temperature. Under the chemical potential difference, water 
molecules and salt ions collide with each other and combine 
together to form small clusters which are the initial nuclei and 
the formation of the nuclei relieves the supersaturation. It is 
difficult to initialize the crystallization in the clear aqueous 
solution; otherwise, the subcooling degree needed for the 
homogeneous nucleation is very large. However, heterogeneous 
nucleation occurs in most cases, in which the dust particles in the 
solution or on the surface of the blade and container will help the 
formation of the nuclei. 

Similar to MPCS, subcooling phenomenon is unwanted due to 
the excess energy consumption in the hydrate generation. 
Several methods have been proposed to suppress the subcooling 
degree, such as adding nucleating agent, agitation, ultrasonic 
vibration and so on, which are similar to the methods used in 
the case of MPCM. Takao et al. [45] suggested to drop fine 
nucleating particles into the solution to act as the nuclei to 
suppress the subcooling, which showed that the subcooling 
degree was reduced from 14 °C to less than 2 °C. Although this 
method can be used to accelerate the nucleation, the nucleating 
agent particles should be prevented from being circulated in the 
entire system, therefore, additional devices are needed to 
reclaim the nucleating agent particles which make the system 
complicated. 

On the other hand, the addition of some nucleating agents, such 
as ethylene glycol and propylene glycol, will change the melting 
temperature of the clathrate hydrate. Wu et al. [46] added 6-8 wt% 
NaCl to 40.5 wt% TBAB solution, which results in decreasing the 
melting temperature from near 12 °C to 6-8 °C. Therefore, this high 
concentration solution can be useful for application due to the 
decrease of the melting temperature. In general, the hydrate 
crystal easily nucleates and grows on the container wall, and then 
the heat transfer performance or energy storage capacity may be 
deteriorated due to the gradual increase of the thickness of the 
hydrate crystal layer. Daitoku and Utaka [47] used the scraper to 
detach the solid layer of TBAB hydrate adhering on the wall which 
would reduce the heat transfer performance, and the different 
morphology of hydrate crystals on the wall before and after 
scraping was reported. 


n<1 
n=1 


Herschel-Buckley 


Bingham plastic n<4 


Newtonian 


Shear stress t 


Dilatant 


Shear rate ý 


Fig. 17. Typical flow curves of purely time-independent non-Newtonian fluids. 


3.3. Flow characteristics of CHS 


3.3.1. Fluidity of CHS 

Different from MPCS, CHS is usually considered as non- 
Newtonian fluid due to its solid—liquid two-phase mixture feature. 
However, it can possibly be treated as Newtonian fluid when the 
solid fraction is very small. 

Darbouret et al. [48] suggested that the traditional viscometers 
were not suitable for measuring the viscosity of the solid—liquid 
suspension due to the density segregation phenomenon. If the 
solid fraction further depends on the temperature, like CHS, the 
measurement is even more difficult due to the temperature 
controlling problem. Consequently, many researchers determined 
the apparent viscosity of suspensions, such as ice slurry, 
microemulsion, CHS, etc., by the shear rate and shear stress which 
are obtained by measuring the volumetric flow rate and 
differential pressure drop of slurry flowing in a horizontal tube. 

Simple non-Newtonian fluid is a time-independent viscous 
fluid, and the shear stress is a function only of shear rate but it is 
more complicated than Newtonian fluid. Fig. 17 shows the flow 
curves of four typical non-Newtonian fluids, including pseudo 
plastic, dilatant, Bingham and Herschel-Buckley fluids. Many 
constitutive equations have been developed to describe the flow of 
non-Newtonian fluid, which are shown as follows 

Power law fluid: 


t=Ky" (11) 


where K is the fluid consistency coefficient and n is called the flow 
behavior index representing the degree that the fluid deviates from 
the Newtonian fluid. n < 1 indicates that the viscosity decreases as 
the shear rate increases, therefore the fluid is shear thinning 
pseudo plastic fluid; n > 1 indicates that the viscosity increases as 
the shear rate increases, therefore the fluid is shear thickening 
dilatant fluid. Power law fluids are with zero yield stress, while 
Eq. (11) is re-written as follows when the yield stress is larger than 
zero, 
Bingham fluid equation: 


T= To + Moy (12) 


Herschel-Buckley fluid equation: 


T = To + Ký” (13) 


The above two equations are the descriptions of Bingham fluid and 
Herschel-Buckley fluid, respectively. Both fluids have yield stress 
To > 0, however, after overcoming the yield stress, Bingham fluid 
behaves as Newtonian fluid and Herschel-Buckley fluid behaves as 
pseudo plastic fluid. 


P. Zhang et al. / Renewable and Sustainable Energy Reviews 14 (2010) 598-614 611 


For the laminar flow in circular pipe, the shear stress varies 
linearly with radius: 
ton e DAP r 

~ YR AR 
where Tw is the shear stress on the pipe wall. Then the volumetric 
flow rate is: 


(14) 


R R 
c= f anrvar =" | rydr= 
0 0 


3 Tw 
= 1 tyýdt (15) 
Ty 0 


The apparent viscosity on the wall is defined as: 


Tw 
Ha(Tw) = Vw (16) 


Shear rate on the pipe wall can be calculated by: 


a E ane Paina E 


din Fea (17) 


"din tw 4 din ty 


where yy Stands for the wall shear rate of Newtonian fluid, 
equaling 4G/7R?, 
Define: 


— dintw 
din Ywn 


Then, for power law fluid: 


__ din(DAP/4l) 
din(8v/D) (18) 


Tw = Ký! =K ý"y (19) 


Therefore, after obtaining the pressure drops and flow rates of 
CHS, the flow behavior index n can be identified by the slope of the 
fitting curve of In(DAP/41) versus In(8v/D). 

For Bingham fluid, the Buckingham equation can be deduced 
from Eq. (15): 


4 
= Tw |1 1 to, 1 /To (20) 
mR Lo 4 3 Tw 12 Tw 
Neglecting higher order item of (to/tw)*, Eq. (20) can be written 
as: 


4 8v 
Tw = 370+ Mop (21) 


Thus, To and 4o can be obtained from the curve of DAP/4I versus 
8v/D. 

Darbouret et al. [48] treated the TBAB CHS as Bingham fluid, and 
applied Eq. (21) to fitting the experimental data. And Hayashi et al. 
[49] measured the flow characteristics of 35, 29, 27 and 26 vol% 
TBAB CHS in the pipes with the inner diameters of 27.6 and 
52.9 mm and the equation for power law fluid was used to fit the 
experimental data. Results of both works are shown in Fig. 18. It 
can be seen from the figure that results of Hayashi et al. seemed to 
vary significantly even with only a small change of volume fraction, 
and the apparent viscosity was larger than that obtained by 
Darbouret et al. at the same volume fractions. 

As shown in Fig. 18, the results of Darbouret et al. indicated that 
there was a critical volume fraction for TBAB CHS at which its flow 
behavior changed from Newtonian fluid to Bingham fluid. 
However, the authors found that it was difficult to determine 
precisely the value of this critical point. Apparently, CHS in higher 
volume fraction has larger yield stress and apparent viscosity. It 
was further found that the viscosity of type B TBAB CHS was a little 
larger than that of type A due to the different shapes of two types of 
hydrate crystals. According to Darbouret et al., the yield stresses of 
both type A and type B CHS varied with ¢°. 

Xiao et al. [50] tested the flow characteristics of TBAB CHS in a 
PVC pipe with 21 mm inner diameter. Power law fluid equation 
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Fig. 18. Flow curves and apparent viscosities obtained by Darbouret et al. [48] and 
Hayashi et al. [49]. 


was applied, which showed that the apparent viscosity of 0- 
16 vol% TBAB CHS was only 3-8 times lager than that of water, 
much smaller than the results obtained by Darbouret et al. [48] and 
Hayashi et al. [49]. And a linear correlation between the flow 
behavior index n and volume fraction @ was proposed based on the 
experimental results: 


n=-2.76+1 (22) 


Flow behavior of CHS may be affected by many factors, such as 
the shape and volume fraction of hydrate crystal particles, the pipe 
shape and size, and so on. It can be concluded from the different 
results given above that the flow characteristics of CHS are very 
complicated and still need more experimental and theoretical 
investigations. 


3.3.2. Friction factor of CHS 

Friction factor can also be calculated through measuring the 
pressure drop and flow rate. For power law fluids, Metzner and 
Reed [51] proposed a modified Reynolds number Rey as follows: 


_ Dvp( 4n \"(8v\'"  D"v?-"p 
Reme =- Ga) (D) Kig" (23) 


Using this modified Reynolds number, the Fanning friction 
factor for the laminar flow of power law fluid can be calculated by 
the following equation which accords with that for the laminar 
flow of Newtonian fluid: 


16 
~ Reyr 


(24) 


for the turbulent flow, only empirical or semi-empirical equations 
were available. Metzner [52] proposed an empirical equation in the 
similar form to Blasius equation for non-Newtonian turbulent 
flow: 


a 
b 
Reir 


f= (25) 


where a and b are the function of the flow behavior index n, and 
both values are given in Table 5. 
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Table 5 

Values of a and b in Eq. (25) [52]. 
N a b 
0.2 0.0646 0.349 
0.3 0.0685 0.325 
0.4 0.0712 0.307 
0.6 0.0740 0.281 
0.8 0.0761 0.263 
1.0 0.0779 0.250 
1.4 0.0804 0.231 
2.0 0.0828 0.213 


Dodge and Metzner [53] also proposed a semi-empirical 
equation for the turbulent flow of power law fluid: 


1 4.0 2-n, 0.4 
JF no7s 08(Reme >) pra oe 


For Bingham fluid in the laminar flow, the following equation 
can be derived from Eq. (20): 


16 1 Het 
Í = Rez (: ) (27) 


3 f?Re3 

where Reg = pvD/[4) stands for the Reynolds number for the 
Bingham fluid, and He = D? pTo/u? is the Hedstrom number. The 
highest order term in this equation can always be neglected. 

Hayashi et al. [49] presented the results of f versus Remr by 
using Eq. (24), and Darbouret et al. [48] gave out the results of AP 
versus Reg based on Eq. (27), and both authors claimed the good 
agreement between the calculated and experimental values. In the 
turbulent flow of CHS, Xiao et al. [50] claimed that both initial mass 
concentration of the aqueous solution and volume fraction of 
hydrate particles had little effect on the results of fand Re, and the 
following empirical equation was obtained from the experimental 
results: 


_1 He 
" 6 Reg 


f =0.3226Re,04271 (28) 


In order to reduce the large frictional pressure drop of CHS flow, 
the drag-reduction surfactants are often added into the fluid. 
Suzuki et al. [54] applied oleylbishydroxyethylmethylammonium 
chloride [C;gH35N(C2H40H)2CH3CCl, trade name: ethoquad 0/12] 
to TME CHS, and found that, however, the agent could only reduce 
the pressure drop of CHS with high solid fraction in the turbulent 
flow, and the heat transfer performance was reduced even more 
because of the suppression of the turbulence. 

There are few results of CHS flow characteristics in the open 
literature. Recent experimental results (not published yet) 
obtained by a group at Shanghai Jiao Tong University (SJTU) 
showed that the flow of TBAB CHS displayed power law fluid 
feature. The flow test was conducted with a horizontal copper tube 
of 6 mm inner diameter. Flow behavior indexes of 5-20 vol% CHS 
were about 0.91-0.74 and varied little when the volume fraction 
further became larger. The experimental results for the turbulent 
flow state were also obtained and showed good agreement with 
Eq. (25) and Eq. (26), as shown in Fig. 19. 


3.4. Heat transfer characteristics of CHS 


Heat transfer tests of TBAB CHS were carried out by the 
researchers in JFE Corporation [55], in which a fan coil unit with a 
nominal flow rate of about 5 kg/min and 3.66 kW cooling capacity 
was used to study the heat transfer characteristics of TBAB CHS. 
The overall cooling load and mass flow rate were measured for 
both chilled water and CHS, consequently, the heat transfer 
coefficient h was obtained. An example of the results is shown in 
Fig. 20. 
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Fig. 19. f-Remr results of TBAB CHS obtained at SJTU. 
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Fig. 20. Heat transfer performance and cooling capacity of CHS in a fan coil heat 
exchanger [55]. 


Compared to the chilled water, CHS exhibited a higher overall 
heat transfer coefficient due to the associated phase change, as 
shown in Fig. 20. However, large viscosity of CHS results in a low 
Reynolds number compared to water at the same flow rate, which 
consequently deteriorates the heat transfer performance. Results 
in Fig. 20 show the cooling capacity of TBAB CHS is about 2.5-3.5 
times larger than that of the chilled water at the flow rate of 1- 
5 kg/min. 

Xiao et al. [56] investigated the heat transfer characteristics of 
TBAB CHS in a horizontal copper tube with 8 mm inner diameter. 
They used the electric resistance as the heater, and the maximal 
heating power was 12 kW. The results indicated that the influence 
of the volume fraction of solid phase on Nusselt number was small 
despite the change of Reynolds numbers. A correlation of Nu with 
Rep was proposed as follows: 


Nu = 0.2049Re0,8732 (29) 


It is obvious that the heat transfer characteristics data is very 
important for the system design and cooling power evaluation of 
CHS, thus, more data is still needed to facilitate this purpose. 

It should be pointed out that the heat transfer of CHS flow in the 
tube is very complicated, which involves not only the variation of 
the volume fraction of solid phase and the consequent variation of 
the melting temperature of CHS but also the change of the flow 
characteristics of CHS because of its non-Newtonian fluid feature. 
The latter two difficulties are not involved in the case of MPCS. 
Furthermore, the fluid flow and heat transfer behaviors of CHS in 
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Fig. 21. The scheme of the air conditioning system using TBAB CHS [57]. 


the heat exchangers obviously urge more investigation on this 
topic as well. 


3.5. Application of CHS to the refrigeration and air conditioning 
system at the industrial scale 


An air conditioning system using TBAB CHS as both the cold 
storage and transportation media in the secondary refrigeration 
loop has been applied to an office building in NKK Keihin 
Workshop [57]. Shown in Fig. 21 is the schematic layout of the 
air conditioning system using TBAB CHS in which an absorption 
refrigerator was used to generate cold. A CHS generator was 
applied to generating TBAB CHS during the off-peak time, and the 
CHS was stored in one of the two tanks, as shown in Fig. 21. TBAB 
CHS was pumped directly to the load side by a special CHS pump 
during daytime to provide cooling, which resulted in the electricity 
peak-shaving. After passing the heat exchanger, TBAB CHS 
dissolved and became the aqueous solution which was then 
stored in another tank. However, if TBAB CHS could not suffice all 
the cold loads, the refrigerator would start to compensate. 

Because of the involvement of CHS in the air conditioning 
system, the initial capital cost was reduced due to the capacity 
reduction of the refrigeration system and the size reduction of the 
pipes, and the operating cost was also reduced due to the reduction 
of the pumping power. The operation during the summer in the 
year of 2000 revealed that the pumping power consumption was 
reduced to 50% of that for chilled water system, and 21% of the 
overall power consumption was reduced. 


4. Conclusions 


An overview of MPCS and CHS used as the phase change 
material slurries (PCS) for the thermal storage and transportation 
has been carried out, which is presented by dividing into two big 
sections. The materials and fabrication methods suitable for MPCS 
and CHS have been presented, and the determination of the 
thermal properties of MPCS and CHS, like thermal conductivity, 
specific heat and viscosity, are also discussed. The flow character- 
istics are different for MPCS and CHS, which can be characterized 
by Newtonian fluid feature for MPCS with small solid volume 
fraction and by non-Newtonian fluid feature for CHS and MPCS 
with high solid volume fraction. The heat transfer of PCS with 
phase change is quite complicated, which still needs more 
intensive investigation. The application of PCS to the secondary 
loop as the thermal storage and transportation media is also 


presented, which is very promising for the practical application in 
the future. 
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